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An optically thin plasma loses more energy to bremsstrahlung than it gains from 

fusion reactions, unless the ion temperature can be elevated above the electron temperature. 
In thermal plasmas, the temperature differences required are possible in small Coulomb 
logarithm regimes, characterized by high density and low temperature. The minimum 
Lawson criterion for thermal plasmas and the minimum pR required for ICF volume 

ignition are calculated. Ignition could be reached more easily if the fusion reactivity can be 
improved with nonthermal ion distributions. To establish an upper bound for this utility, 
we consider a monoenergetic beam with particle energy selected to maximize the beam- 
thermal reactivity. Channeling fusion alpha energy to maintain such a beam facilitates 
ignition at lower densities and pR, improves reactivity at constant pressure, and could 
be used to remove helium ash. The gains realized with a beam thus establish an upper 
bound for the reductions in ignition threshold that can be realized with any nonthermal 
distribution; these are evaluated for p-^^B and DT plasmas. 


I. INTRODUCTION 

Fusion reactions which release most of their 
energy in charged particles are desirable for 
power applications. In particular, if the num¬ 
ber of energetic neutrons producecPis small, the 
power plant can be designed with less expensive 
shielding and with fewer material constraints. 
Also, direct conversiorP of charged particle en¬ 
ergy could offer such a scheme much greater ef¬ 
ficiency than a thermal cycle. The p-^^B reac¬ 
tion is ideal due to the low incidence of neutron¬ 
generating side reactions and the reactants’ nat¬ 
ural abundance. However, at typical densities, 
an optically thin p-^^B plasma loses more en¬ 
ergy via bremsstrahlung radiation than it gains 
from fu sion reactions, making sustained burn 
difficult .1^ Lacking sustained burn, or what is 
referred to as ignition, does not mean that fu¬ 
sion energy cannot be extracted, as in a wet 
wood burner. But the lack of ignition makes 
such a means of extracting energy technologi¬ 
cally difficult and expensive.^ 


Electronic mail: hay@princeton.edu 


Previous efforts toward p-^^B ignition have 
focused on mitigating the bremsstrahlung loss 
channel. One approach is embedding large 
magnetic fields in the fusing plasma, restrict¬ 
ing the motion of unbound electrons and re¬ 
ducing the bremsstrahlung emission both di¬ 
rectly, by modifying electron-electron scatter¬ 
ing, and indirectly, by modifying electron-ion 
scattering. The Landau wavefunctions differ 
considerably from the isotropic, field-free case 
and the electron-electron scattering assumes a 
I-D character. In the aggregate, the spectrally 
integrated emission is reduced by about 20%.!^ 
Likewise, by increasing the Landau energy level 
spacing {Ml), a hot ion mode {Ti > T^) can 
be preserved due to the suppression of ion- 
electron energy transfer^ such a mode is char¬ 
acterized by low bremsstrahlung emission at a 
given plasma pressure. 

In a plasma with substantial electron degen¬ 
eracy, bremsstrahlung is also reduced. The 
main effect appears to be a reduction of the 
electron stopping power,^^ enabling a large 
electron-ion temperature difference. For the 
same reason, degenerate electro ns are prefer¬ 
able in beam-initiated detonation^^IIoIII] of com- 
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pressed fuel that rely on rapidly heating 

the ion species in a small locus. 

Ignition using is possible when the ion 

temperature is raised significantly above the 
electron temperature, a circumstance realized 
only in high-density thermal plasmas. (Though 
this does occur naturally in low temperature 
plasmas because of electron degeneracy.) Ti 
Tg is necessary, but not sufficient: Ti itself must 
be on the order of several hundred keV due 
to the small p-^^B fusion cross section below 
that threshold. The resulting plasma is strongly 
collisional, but the robust fusion burn sustains 
ion temperatures above the electrons cooled by 
bremsstrahlung (note that \Ti — Tel ^ Pfl^ie 
under these conditions). 

Because high densities and temperatures are 
difficult to reach simultaneously, it is important 
to consider how either of these requirements 
might be relaxed. In equilibrium ignition- 
regime p-^^B plasmas, fusion alpha particles 
slow chiefly on the ions. Using waves to channel 
alpha energ}112*i^ to the lighter fusing species 
(protons) could improve the reactivity at fixed 
pressure, but, in the case of p-^^B as opposed to 
DT, the net gain is limited by the small fraction 
diverted to electrons in equilibrium. 

Nonthermal distributions offer a plausible 
means to ease ignition requirements, either 
by im provin g the MHD stability of the bulk 
plasma,fi^E^ or by increasing the number of re¬ 
acting particles near the cross section peak in 
velocity space, as treated here. This work con¬ 
siders the effect of such a non-Maxwellian fea¬ 
ture in the ion distribution, namely a monoen- 
ergetic beam with particle energy maximizing 
the beam-thermal reactivity. 

In order to discover the minimal ignition cri¬ 
teria, it is assumed that a fraction of the fusion 
power is efficiently channeled to maintain the 
beam against collisions with thermal particles. 
By way of comparison, these ignition conditions 
are calculated for both DT and p-^^B plasmas, 
with and without the presence of the monoener- 
getic beam, in both magnetic and inertial fusion 
configurations. 

In all instances, power flows from a ‘hot’ 
species (alphas, protons) to a cool one (pro¬ 


tons, boron, respectively), and the extra power 
needed to maintain the nonthermal distribu¬ 
tions is consistently included in the calculation. 
We explicitly include the effects of ion-ion and 
electron-ion drag in the energy balance. 

However, we do neglect the power flow re¬ 
quired to maintain the beam against velocity 
space diffusion. Consider a test particle born 
with energy slightly larger than that associated 
with the fusion cross section resonance. In both 
p-^^B and DT plasmas, this particle’s fusion 
rate competes chiefly with slowing on light ions 
and electrons (boron excluded). Only after the 
particle slows out of the fusion resonance do the 
collisional drag and velocity space diffusion due 
to thermal ions become significant. Because the 
diffusion is especially small compared to drag at 
the beams’ typical energies, we are justified in 
neglecting this term in the power balance (see 
Appendix A). 

Likewise, we do neglect kinetic effects such 
as collisionless streaming instabilities. A cold 
beam in a warm plasma will be subject to a 
complex array of instability mechanisms, but 
the two-stream instability between the fast mo- 
noenergetic beam ions and warm bulk species 
is prognostic.^ However, in both DT and p- 
^^B ignition-regime plasmas, the bulk species 
are sufficiently cool with respect to the beam 
that nearly monoenergetic beams are resilient 
to this instability (see Appendix 

Likewise, favorable assumptions and esti¬ 
mates are invoked to position the calculation as 
an upper bound on the feasibility of these fusion 
scenarios. For example, we explicitly assume 
that free energy in the alpha distribution can 
be transferred to the protons with perfect effi¬ 
ciency, whereas in practice rf power will likely be 
required to establish and maintain the mediat¬ 
ing waves against collisions, etc. However, these 
beams do provide an upper bound on the gross 
gains (cf. net gains which would include the cost 
of maintaining a self-consistent distribution) re¬ 
alizable from nonthermal distributions; each ad¬ 
ditional particle in the beam adds the maximum 
amount of reactivity that could be gained from 
an extra particle at any energy. 

In all cases, the relative locations of the ther- 
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mal and beam reactivities’ maxima determine 
the possible advantages of a nonthermal scheme. 
In particular, a successful beam contributes ex¬ 
cess reactivity at constant pressure, so a plasma 
burning at a temperature close to the peak cross 
section energy (as in the case of DT) reaps 
limited benefit from nonthermal features. Al¬ 
though the igniting operating regime is difficult 
to access, the maxima of the thermal and 
beam reactivities are germane to a nonthermal 
configuration which substantially lowers the ig¬ 
nition threshold. This work suggests the exis¬ 
tence of nonthermal regimes where igni¬ 

tion may be possible. 

In section II, the model used to analyze the 
equilibrium state of a nonthermal plasma is de¬ 
tailed. Sections III and IV apply this model 
to predict ignition thresholds in magnetically- 
and inertially-confined plasmas. Section V con¬ 
siders the potential benefits of non-Maxwellian 
features in these plasmas. 

II. THERMAL EQUILIBRATION MODEL 

In the case of p-^^B , the proton population is 
modeled as a thermal bulk plus a fast monoener- 
getic beam located near the p-^^B fusion cross 
section peak at 592 keV.^^^ A 0-D equilibrium 
model provides self-consistent species tempera¬ 
tures: 


PbPfus = PsD 

(la) 

PePfus + Vgp{Tp — Tg) 
+VgB{TB —Te)— Phrem = 0 

(lb) 

PpPfus + BpgiTg — Tp) 
+Vpb{Tb — Tp) = 0 

(Ic) 

PBPfus + i^Bp{Tp — Tb) 

(Id) 

+VBe{Tg — Tb) = 0 


pPi Bb^Bp^TlB^ Tg^Tp^TB) 

(le) 


Eqs. (la-d) describe the flow of fusion charged 
particle power to each species; Eq. (le) is in¬ 
dependent of the first four and determines the 
volume gain of an assembly with temperatures 


and densities characterized by (la-d) once the 
scale of the system (p) is specified. 

An effective frequency Vij = Vji describes the 
drag of a thermal population i on thermal pop¬ 
ulation These coefficients are notably sen¬ 
sitive to the Coulomb logarithms characterizing 
the interacting species; the plasma electron den¬ 
sity is an important parameter in establishing 
ignition criteria. The fusion power is the sum 
of a thermal reactivitjC^ and a beam-thermal 
reactivity,^ viz. Pf^^/Wf = nsrip {crv)^^ + 
riBrib {crv)i^^, where Wj = 8.7MeV is the en¬ 
ergy released in one reaction. The thermal 
bremsstrahlung emission Pbrem is calculated as 
that from an optically thin medium, including 
relativistic corrections up to {Tf,/meC^yW^ 

The Pi denote the fraction of the fusion power 
Pfus deposited in the species, such that 
Pi = 1. In particular, Pb is the fraction 
of the fusion power spent preserving the pro¬ 
ton beam velocity distribution, such that Pb = 
Psu/Pfus, with PsD = Y.i'^^bVbiEb the amount 
of power dissipated by fast protons slowing 
down in the plasma (see Appendix A). The 
Vbi ’s are effective frequencies describing the fast 
beam particles’ energy loss to the i**' species.l^^ 
The constant beam particle energy Eb is dis¬ 
placed from the cross section peak to maximize 
the beam-thermal reactivity at a specified Tg. 
In order to limit the anticipated damping of 
waves used to construct the beam, Eb was re¬ 
stricted to values at least twice as large as Tb- 
In the large Tg limit anticipated here, fusion 
alpha particles slow chiefly on the ions, and the 
amount of alpha power diverted to electrons is 
of order 10%. In order to estimate the amount 
of alpha power absorbed by each species, con¬ 
sider a fast alpha particle born in a thermal p- 
^^B plasma with energy Ei. The amount of that 
particle’s energy deposited in the protons is 

7’p ~ J^apE{t) dt (2) 

where E{t) is the instantaneous alpha particle 
kinetic energy, v^i is the effective energy loss 
frequency for fast alpha particles slowing on the 
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species. Ef = E[tf) is the largest alpha en¬ 
ergy at which energy flow from the particle to 
one species of field ions equals that from the 
field ions to the particle. As the typical alpha 
particle slows down from « 3MeV in a 
plasma with Tp = Tb = it reaches this dy¬ 
namic equilibrium with the electrons first: al¬ 
though net energy exchange with the electrons 
has ended, the alpha particle continues to heat 
the proton and boron distributions at the tf de¬ 
fined this way. Thus, Ef Ki 1.5Te, such that 
Vae{tf) vanishes. Changing the variable of inte¬ 
gration to the particle energy, note 

(3) 

i 

and find 

Ei 

J ^ap “T ^aB “r ^ae 

Ef 

ifB and ipe are defined analogously. This model 
satishes (fiB + <Pp + P’e = 1; the distribution of 
fusion product energy is sensitive to the three 
species’ temperatures as well as their relative 
concentration, determined by e = ua/bb, the 
number ratio of the heavy species ‘A’ to the 
light species ‘B’ (here e = n^/np). 

In the model problem of an equilibrium p-^^B 
reactor, where a fraction Pb of the fusion power 
is devoted to counteracting the proton beam 
slowing down, only a fraction /3p = {1 — /3b)^p 
is available to the protons. This /3b is extracted 
from alpha particles near their birth energies of 
several MeV. Because a temperature-dependent 
spectrum of birth energies peaking around 3.5 
MeV is observed for the p-^^B reaction, it is as¬ 
sumed Ei = 3 MeV qualitatively predicts the 
energy flow to the thermal ions.^^ The model 
includes species absorption fractions calculated 
accordingly. 

III. MAGNETIC CONFINEMENT 


species temperatures and the amount of power 
spent preserving the monoenergetic proton 
beam. Because the fusion reaction rates and 
bremsstrahlung emission both depend on 
only the ratios of species densities (rib/np and 
e = ubIbp) are important. However, the tem¬ 
perature relaxation rates are sensitive to the 
Coulomb logarithm log A, which in turn de¬ 
pends on the absolute electron density. 

In order to explore the full space of igniting 
plasmas, consider a wide range of electron den¬ 
sities (10^° < He < 10^°cm“^) and disregard 
the boron temperature equation (Id), instead 
taking Tb as a specified parameter. The set 
(la-c) can be solved for /3b, T/., and Tp once rie, 
Tb, nb/up, and e are chosen. This specihcation 
of Tb is tantamount to asserting an arbitrary 
ion energy confinement time. 

With these assumptions, ignition thresholds 
can be defined for magnetically conhned plas¬ 
mas. Lawson criteria are presented for thermal 
p-^^B and DT plasmas, followed by nonthermal 
p-^^B and DT plasmas. These latter calcula¬ 
tions provide the upper bounds on the utility of 
maintaining a nonthermal distribution, abiding 
by the assumptions laid out in the Introduction. 
In the thermal cases, the numerical model pre¬ 
dictions are compared to analytic calculations. 

A. Lawson criterion 

A common metric for the performance of fu¬ 
sion power systems is the criterion first ob¬ 
tained by Lawson.l^ In steady state, a plasma 
ignites if the nr product of number density and 
conhnement time exceeds a specific value de¬ 
pending on the fusion reaction under consider¬ 
ation. The Lawson criterion is the statement 
of this minimum value, determined simply by 
Pius > P\oss- Taking Pfus = WfUAUB 
with Wf the energy released in charged parti¬ 
cles, and defining r = Wj{—dWjdt) where W 
is the thermal energy content of the plasma and 
-dW/dt = Pbrem = Ploss^ 


The four equilibrium equations (la-d) de- rieT > + Zb + e(l -I- Za)) 

scribe the time rate of change of the ^ 
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where e = nA/ns, ne = ZAnA+Zuns, and T is 
the temperature common to all species. In gen¬ 
eral, the species temperatures differ, rendering 
the minimization of the right hand side nontriv¬ 
ial. It is however instructive to search numeri¬ 
cally for those operating conditions which afford 
the least stringent Lawson criterion. The re¬ 
sulting numerical criterion is necessary but not 
sufficient for ignition. In pressure-limited sys¬ 
tems, the triple product nTr is a superior met¬ 
ric because it is proportional T^/ {ov), in turn 
inversely proportional to the achievable fusion 
power W fp^ {(Tv) /T^.I^Thus the threshold ig¬ 
niting state has minimum nTr and maximum 
power. 

B. Minimum Lawson criterion for DT ignition 


the extremal point at a mass-weighted ion tem¬ 
perature Tion = 8.7 keV and e = 0.97. The 
corresponding minimum value of nTr is 2.6 x 
10^^ keV• s/m^, in close agreement with the es¬ 
timated value. A slight excess of deuterons is 
explained by their more efficient use of thermal 
energy. At fixed energy, the lighter particles 
have larger velocities, and for a thermal 9 keV 
DT plasma, larger cross sections (cr/ is mono¬ 
tone increasing up to Ecm = 64keV). Balanc¬ 
ing this effect are the fusion rate penalty with 

Plus OC (8) 

and deuterons’ larger drag losses to the elec¬ 
trons; the ratio of temperature equilibration 
rates is voe/^Te ~ 1-5. 


Assuming Tg = T^ = T, one can make a 
simple estimate of the threshold Lawson crite¬ 
rion for DT ignition. Because both deuterium 
and tritium carry the same number of electrons, 
Tbrem docs not depend directly on the number 
ratio e. The fusion power, however, is maxi¬ 
mized for e = nx/nD = 1, so the extremal ig¬ 
nition case will be located near that operating 
point. Using Eq. (§, we can form the triple 
product criterion for DT: 


nTr > 


2T^ (1 -H ef 

Wf {(Tv) e 


/(e)5(r). 


( 6 ) 


In other words, the minimum triple product is a 
separable product of functions of e and T. Con¬ 
sidering minimization with respect to e, we have 


min(nTT) (x 2 -|- e H—, 


(7) 


which has its minimum at e = 1. The crite¬ 
rion is also proportional to T'^ f{av), a function 
of temperature only. The numerical optimum 
temperature should be located close to the cri¬ 
terion’s minimum at T = 13.5 keV. The esti¬ 
mated lower bound is 1.9 x 10^^ keV- s/m^. 

Alternately, the consistent set given by 
Eqs. (la-d) can be solved at various Ti to de¬ 
termine the threshold (minimum nTr) ignit¬ 
ing configuration. A numerical search found 


C. Minimum Lawson criterion forp-^^B ignition 

Eollowing the approach of the previous sec¬ 
tion, the triple product criterion has the density 
scaling 

min(nTr) oc 8 H-h 15e, (9) 

e 

which takes its minimum value at e = ns/np = 
l/-\/l5 « 0.26. Likewise, the quantity j{av) 
is minimized for T = 138 keV. This inconsistent 
optimization, regarding e and T as totally inde¬ 
pendent quantities, returns an estimated lower 
bound of 5.1 x lO^^keV- s/m^. However, a p- 
^^B plasma cannot ignite with = Ti. (That 
is, unless a very substantial number of x-rays 
can be reflected off the walls and reabsorbed by 
the plasma. Here, we assume that for all prac¬ 
tical purposes, this cannot be done. We like¬ 
wise assume that the plasma is optically thin 
to bremsstrahlung.) Per the analytic model, 
Tfus/Tbrem kas a maximum value of 0.44 at 
Tg = T, = 204 keV and e = 0.11. 

A numerical search allowing for species- 
dependent temperatures found the extremal ig¬ 
nition point at a mass-weighted ion temperature 
Tion = 193 keV and e = 0.26. The correspond¬ 
ing minimum value of nTr is 2.2 x lO^^keV- 
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s/m^, in close agreement with the estimated 
value. The numerical optimum is cooler and sig¬ 
nificantly more boron-rich than the naive min¬ 
imization of Eq. (§ would suggest. The ex¬ 
tremal ignition state offsets a larger boron con¬ 
centration (more fusions and bremsstrahlung) 
with a cooler = 59keV. 


D. Nonthermal gains: DT 

Suppose there were some way, say by alpha 
channeling, to support a monoenergetic deu¬ 
terium beam. In both the DT and plas¬ 

mas, the lighter species is chosen for the beam 
because of the lower energy investment required 
to achieve the high center of mass energy needed 
to access the fusion cross section peak. Oper¬ 
ating points with A = log A < 3 were discarded 
because the coupling coefficients are accu¬ 
rate only to first order in A“^. Points with 
/3f, > 1, signaling the need for injected power 
to maintain the beam, were discarded likewise. 

The maximum reactivity subject to a pres¬ 
sure constraint is a good metric for comparing 
the effects of nonthermal distributions in fusion 
plasmas because the pressure a reactor can con¬ 
fine is limited by magnet strength. In the case 
of ITER, this figures to be about 10 bar, which 
we adopt as a standard value in order to com¬ 
pare DT and MFE schemes. Nonther¬ 

mal features which improve the fusion reactiv¬ 
ity at a fixed or reduced pressure are therefore 
desirable. In general, the plasma pressure in¬ 
cludes contributions from any beam as well as 
the thermal populations and any alpha parti¬ 
cles slowing down on thermal particles. Because 
the alpha particles whose energy is channeled to 
maintain the beam (fraction are ‘lost’ on a 
fast, collisionless timescale, their pressure can 
be assumed to be negligible. The total alpha 
pressure is estimated as 

(l-/3,K^ = (l-/3,)-^ (10) 

Na Nal^SD 

where Ua is the number density of alpha par¬ 
ticles, No, is the number of alpha particles 
spawned by a single reaction, and vsb is the 


slowing down frequency on thermal particles. 
In practice, wave-mediated diffusion and device 
confinement could supersede Coulomb collisions 
with thermals as the salient processes limiting 
the average alpha lifetime. The reactivity in¬ 
cludes contributions from both thermal-thermal 
reactions and beam-thermal reactions. 

The ion number ratio which maximizes 
the reactivity of a pressure-limited, constant- 
temperature system (in the absence of a beam) 
is (see Appendix B) 


^0 — + Zb )/{1 + Za )- ( 11 ) 


The temperature Tq at the reactivity maximum 
satisfies 


A 


{av){T) 


T=To 


^ {av){To) 

To 


( 12 ) 


Note that both of these conditions are indepen¬ 
dent of the limiting pressure. In the case of 
a constant-temperature DT plasma, the opti¬ 
mal reactivity is found at T = 13.5 keV and 
e = 1. In the limit of small beam fraction 
(fi —>■ 0, the numerical model locates the reactiv¬ 
ity optimum under a 10-bar pressure constraint 
at Ti = 15.3 keV and e = 0.89. 

The introduction of a fast deuterium beam at 
the cross-section peak can improve the fusion 
reactivity at constant pressure. Fig. 1 indicates 
the possible gains up to a beam fraction of 1/2. 
The system pressure is held fixed at 10 bar. 

As the beam fraction increases, the 
reactivity-optimizing temperature decreases 
and the ion ratio tilts toward the heavier target 
species. Both of these shifts tend to increase 
the beam reactivity at the expense of thermal 
reactivity. Such an optimized beam reduces 
the necessary thermal energy content of an 
igniting system because fewer fast particles are 
off-resonance; at constant pressure, a larger 
reactivity is possible. Likewise, a fast beam 
reduces the density and pressure required for 
an igniting plasma, substantially easing the 
minimum ignition conditions. 





7 




FIG. 1. Optimized DT operating points with in¬ 
creasing nonthermal features in the light ion distri¬ 
bution function. The deuterium beam fraction 
varies from 10~^ to 0.5; the optima at larger beam 
fractions are characterized by greater reactivity (a) 
[cm®/s], increased target ion concentration (b), and 
lower ion temperatures (c) [eV]. 


E. Nonthermal gains: 

In the beam-free, Maxwellian case, subject 
to a pressure constraint, the reactivity 

is maximized at T = 138 keV and e = 1/3. 
In the limit of small beam fraction, the reac¬ 
tivity maximum (under a 10-bar pressure con¬ 
straint) occurs at = 137 keV and e = 0.20. 
The discrepancy in e is due to the signifi¬ 
cant bremsstrahlung emission present in the full 
model at higher values of e. The possible gains 
from the addition of a fast proton beam are in¬ 
dicated in Fig. 2, up to a beam fraction of 0.1. 

The energy cost of maintaining a given beam 
fraction against collisions with thermals is com¬ 
paratively greater in than DT, due pri¬ 

marily the abundance of electrons. Above about 
ip = 0.2, a proton beam cannot be maintained 
without injected power. In fact, above beam 
fractions of 0.02, these beams were only useful 


FIG. 2. Optimized operating points with in¬ 

creasing nonthermal features in the light ion dis¬ 
tribution function. The proton beam fraction p 
varies from 10“^ to 10“^; the optima at larger beam 
fractions are characterized by greater reactivity (a) 
[cm®/s], increased target ion concentration (b), and 
higher ion temperatures (c) [eV]. 


at higher temperatures, where the slowing down 
of beam protons is reduced, so the beams can be 
maintained at lower cost. Likewise, in the ab¬ 
sence of a pressure constraint, fast beams can 
be used quite profitably. At sufficiently high 
temperature and density, the power required to 
counteract beam drag is a small fraction of the 
total thermal power. 

The removal of fast alpha particles in the 
channeling process is a significant boon to p-^^B 
ignition prospects. In a pressure-limited device, 
the enormous number of fusion alpha particles 
produced poisons the reaction by contributing 
pressure without also providing reactivity. Even 
a marginal increase in the reactant pressure, as 
would be the case with a small beam fraction 
of fast particles, could be helpful because the 
fusion power scales as p^. 

Consider the timescale defined by T = 
Pth/ Pa, where pth is the total pressure of all 





























thermal particles in the plasma and Pa = iana 
is the instantaneous fusion power released in al¬ 
pha particles. If the pressure pth varies slowly 
on the timescale T, for example in the case of 
continuous refueling, T is a reasonable estimate 
of the time a plasma can burn in quasi-steady 
state conditions (for which alpha poisoning is 
insignificant). The ratio 7dt/7),b is of interest. 
Assuming equimolar DT with Tp = It = Ti, 

udTd + rirTr + UeTe 

/DT = -—r- 

CaUDriT [crv) 

^ 4(r, + Te) 
eaUe {av) 

Supplying projected ITER plasma parameters, 
Tdt ~ 1 s, which indicates that controlling al¬ 
pha pressure will be important during pulses 
expected to last several minutes. In compari¬ 
son. 


(13) 

(14) 


TpB = 


Uplp ~\~ TIbTb ~\~ TleT^ 

SeaUpUB (crv) 


= (l + 5e) 


(Tp + Tb + (1 + 5e)Te) 
SCaUe (ctu) 


(15) 


(16) 


where Scq, = 
11 


8.7 MeV. Note that in the case 
of the reaction. Pa = Pins- Supply¬ 

ing the parameters minimizing the Lawson cri- 
terio n for a thermal plasma (see sec¬ 

tion IIIC), TpB ~ 10“^® s. Clearly an active 
means of removing alpha pressure (on a colli¬ 
sionless timescale) is crucial for any plausible 
p-^^B reactor. 

Although the potential gains of a fast beam 
in p-^^B plasmas (due to the high reactivity of 
a beam near the thermal bulk) are limited by 
drag on the electron densities required, alpha 
channeling may yet prove an invaluable means 
of controlling the alpha poisoning effect. 


IV. INERTIAL CONFINEMENT 

Here we describe the volume ignition 
schemd^ of inertial fusion and discuss how DT 
and p-^^B plasma parameters might be opti¬ 
mized to lower the ignition threshold. Although 


conservative in its predictions (cf. practical ICE 
schemes), the volume model is useful because of 
its simple structure (nearly neglecting hydrody¬ 
namic motion), highlighting the effect of a fast 
beam on ignition conditions. In particular, we 
will determine to what extent fast beams could 
reduce the assembly energies and pR of DT and 
p-^^B volume targets. 

The gain equation (le) is independent of the 
equilibrium equations (la-d). To solve the gain 
equation, however, the absolute density p and 
scale (pi?, or equivalently Eg) of the system 
must be specified. Together with the tempera¬ 
tures and beam density supplied by Eqs. (la-d), 
these are sufficient to evaluate the volume gain 
(the yield of a homogeneous spherical assembly 
burning in a sound time divided by its initial 
thermal energy). 


A. Volume ignition 


The volume ignition schem^^ imagines a 
spherical target that has been prepared in a 
completely homogeneous state at the time of ig¬ 
nition. This provides a conservative estimate of 
the obtainable gain because the entire fuel must 
be heated; more tractable ICE schemes rely on 
the heating of only a small portion of the burn¬ 
ing mass, reducing the total thermal energy of 
the assembly, Eg- In practice, the driver energy 
required to assemble the target is greater than 
Eg due to backscatter, x-ray conversion losses, 
solid angle effects, rocket efhciency, etc. 

The gain associated with a volume-ignited 
target can be expressed simply as (cf. Eq. (le)) 


G = 


Eif) 

frTeff 


(17) 


where E is the fusion energy released per unit 
mass (3.39 • 10“ J/g for DT, 7 • 10i° J/g for p- 
^^B ) and (p is the fraction of the target mass 
that burns before hydrodynamic disassembly, 
r = ^B/w is the plasma’s specific gas constant 
(m being the mean mass of the constituent par¬ 
ticles) and Teff is the number-weighted temper¬ 
ature. frTeff = Eg/Mf, where Mf is the total 
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target mass. It is assumed that classical statis¬ 
tics suffice to describe each of the target’s con¬ 
stituent species. 

In order to operate in a high-gain regime, a 
typical IGF target will burn a significant frac¬ 
tion of its fuel.l^ It is therefore necessary to in¬ 
tegrate the fusion rate equation h ~ over 
the confinement time to estimate the fraction 
of fuel consumed in fusion reactions. The result 
for a mixture of two species with initial number 
densities hao and ubo, such that e = nAo/nso, 
is (see Appendix C) 


termined by Eqs. (la-d): 

= = ( 20 ) 

3 Mj[T, + Tp + e{5T,+TB)] 

2 rUe + mp + e{5me + niB) 


in the case of p-^^B with UBo/n-po = e and 
i^b/np = 0 (no beam). Mf is the mass of the en¬ 
tire target; Boltzmann’s constant kB has been 
suppressed. Using 




( 22 ) 


e-f 1 


1 - 


e- 1 


e+l Mb') ^ 


(18) 


where Hb = 6csm/{av), Cg = is the 

ion sound speed, and m = {erriA + rnB)/{^ + 1). 
p = mAnAo + iTT-BnBo is the initial mass density. 
In the limit e —> 1, the familiar pR formula for 
an equimolar target is recovered: 


pR 


pR Hb 


(19) 


the burnup (j) = 4>{pR,Te,Tp,TB) and there¬ 
fore the volume gain are entirely specified by 
Ea once the assembled density p has been given 
(in addition to e and /3b). To wit. 


pR 


1 / n 

/ p^Ea me + mp + e(5TOe -I- UIb) \ 

Te + Tp + e(5Te + Tb) J 

(23) 


The addition of a fast proton or deuteron beam 
is reflected by 

Ea ^ ^ (jlb^b + NgTg^ (24) 


In the limit pR —> oo, burnup is limited to 

2/(e + !)■ 


B. Optimization procedure 

Along with e, the initial ratio of boron to hy¬ 
drogen nuclei, the proton beam fraction rib/up 
is regarded as a parameter. With e and rib/rip 
specified, the system of Eqs. Q can be solved 
for Te, Tb, Tp, and rib- These state variables are 
used to evaluate the volume gain for the config¬ 
uration, which is then a function of pR only 
(through the burn fraction, (/>). A large sam¬ 
pling space was considered to locate the mini¬ 
mum assembly energy Ea necessary for a vol¬ 
ume gain of unity (Eq. (Hzl))- Ea can be re¬ 
lated directly to the equilibrium conditions de- 


Thus, once equations (la-d) have been solved 
consistently for the species temperatures with 
specified Tb, e, beam fraction, and density, the 
gain equation (le) can be solved for the mini¬ 
mum assembly energy resulting in a volume gain 
of unity. Using the above definitions, it is pos¬ 
sible to express this criterion in terms of the 
assembly pR. 

C. DT ignition criterion 

Assuming a constant temperature T for all 
species, it is straightforward to estimate a 
beam-free best case from the gain equation (le). 
At p = 10^, the result is a minimum ignit¬ 
ing pR of 0.022 g/cm^ at T = 18keV and 
e = tit/ no = 1.10. This is equivalent to an 
assembly energy of 0.09 J. 
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A numerical search allowing for distinct, 
consistent species temperatures determined by 
the system (la-e) located a minimum pR of 
0.07 g/cm^ at a mass-weighted ion temperature 
of 7.1 keV and e = 1.13. This is equivalent to 
an assembly energy of 1.1 J. In the full model, 
slowing of fusion alpha particles on the ther¬ 
mal electrons creates some separation in elec¬ 
tron and ion temperatures (Tg « lOkeV here), 
rendering the assumptions of the original esti¬ 
mate inaccurate. 

In the presence of a fast deuteron beam, the 
ignition criteria are further relaxed. The typ¬ 
ical minimal state features a fast beam collid¬ 
ing with a thermal plasma substantially cooler 
than the Maxwellian optimum. In DT, the po¬ 
tential utility of the beam is limited by the fact 
that the resonant beam energy is an order of 
magnitude greater than the optimal Maxwellian 
plasma temperature whereas the reactivity gain 
is limited to a factor of about 3.5. As Fig. 3 
indicates, a fast deuteron beam nonetheless re¬ 
duces the pR required for ignition by a factor 
of 14, corresponding to a factor of 10"^ reduc¬ 
tion in the total thermal energy of an igniting 
assembly. Above a beam fraction of about 4%, 
the shift towards large tritium concentrations 
(Fig. 3b) and low ion temperatures (Fig. 3c) is 
pronounced. 

D. ignition criterion 

Assuming a constant temperature T for all 
species, it is straightforward to estimate a 
beam-free best case from the gain equation (le). 
At p = 10^, the result is a minimum igniting pR 
of 9.95 g/cm^ at T = 153 keV and e = nslrip = 
1.30. This is equivalent to an assembly energy 
of 59 MJ. 

A numerical search allowing for distinct, 
consistent species temperatures determined by 
the system (la-d) located a minimum pR of 
2.1g/cm^ at a mass-weighted ion temperature 
of 196 keV and e = 1.8. This is equivalent 
to an assembly energy of 0.28 MJ, showing 
marked improvement over the naive estimate. 
At the high densities characteristic of ICF, the 



FIG. 3. Optimized DT operating points with in¬ 
creasing nonthermal features in the light ion dis¬ 
tribution function. The deuteron beam fraction (p 
varies from 10“^ to 0.5; the optima at larger beam 
fractions are characterized by lower minimum ignit¬ 
ing pR (a) [g • cm“^], increased target ion concen¬ 
tration (b), and lower ion temperatures (c) [eV]. 


Coulomb logarithm is small enough that a large 
temperature difference can be sustained be¬ 
tween electrons and ions. In the particular case 
of , the substantial bremsstrahlung emis¬ 

sion keeps the electron temperature low. The 
resulting lower thermal content of the fusing 
plasma improves the volume gain. 

In the presence of a fast proton beam, the 
ignition criterion is relaxed. In contrast to 
DT, the resonant beam energy (c. 592 keV) 
is only a factor of 2-4 greater than the opti¬ 
mal Maxwellian plasma temperature, whereas 
the reactivity gain is superior, about a factor of 
6. Fig. 4 suggests that in , the net ef¬ 

fect of a beam is nevertheless less pronounced 
cf. DT. The minimum igniting pR is reduced 
by nearly a factor of 3 (Fig. 4a), corresponding 
to a factor of 30 reduction in assembly energy. 
The optima eschew contributions from thermal 
reactions above a beam fraction of about 7%, as 
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FIG. 5. R metric plotted for constant-temperature 
DT (solid line) and (dotted) plasmas. The 

plotting intervals are restricted by the validity of 
thermal reactivity fits and the incursion of the res¬ 
onant beam into the bulk of the distribution func¬ 
tion. We allow Cbeam > ST. 


FIG. 4. Optimized operating points with in¬ 

creasing nonthermal features in the proton distri¬ 
bution function. The proton beam fraction ip varies 
from 10“^ to 0.5; the optima at larger beam frac¬ 
tions are characterized by lower minimum igniting 
pR (a) [g • cm“^], roughly equal target ion concen¬ 
tration (b), and lower ion temperatures (c) [eV]. 


demonstrated by the shift to larger e (Fig. 4b) 
and low ion temperature (Fig. 4c). 


V. DISCUSSION 

This work seeks to address, in broad strokes, 
the utility of non-Maxwellian features in both 
DT and plasmas under magnetic and in¬ 

ertial confinement. Regardless of the scheme 
chosen, the relation between the thermal reac¬ 
tivity and resonant features in the fusion cross 
section is crucial in establishing this utility. 

A natural criterion in both confinement 
schemes is the amount of pressure or thermal 
energy required for ignition because this quan¬ 
tity scales directly with the facility cost (mag¬ 
net strength in MCF, driver energy in IGF). If 
a beam is to reduce the pressure requirement, it 
should provide excess reactivity without a dis¬ 


proportionate contribution to the system pres¬ 
sure. In particular, the beam-thermal reactivity 
ratio and the ratio of beam energy to the bulk 
temperature should be compared. Heuristically, 
the parameter 


R = 


(^^)beam / 

ebeani/F 


(25) 


is indicative of the beam utility for specified 
plasma conditions. 

In the case of DT, the fusion resonance is 
located far in the tail of most igniting plas¬ 
mas (true ‘thermonuclear’ fusion). Because 
bremsstrahlung losses are minimal, the plasma 
energy content required for ignition is small 
(with respect to the resonance), but the dis¬ 
tance of the resonance is a drawback when en¬ 
ergy is supplied directly to heat particles at 
such large energies, as in the case of a reso¬ 
nant beam. Likewise, the broad DT resonance 
limits the achievable reactivity gain; the value 
of R in igniting thermal DT plasmas is about 
1/2, indicating that channeling fusion power to 
a beam is not generally useful. Fig. 5 traces 
R for constant-temperature DT and p-^^B plas¬ 
mas. 

In contrast, the p-^^B cross section energy 
peak is close to the temperatures relevant to 
thermonuclear fusion. Likewise, due to the nar- 
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row resonance, substantial reactivity is gained 
by channeling fusion power to a resonant beam 
(e.g. the beam reactivity is 15.7X greater than 
thermal at T = 100 keV). Typical values of R 
in an igniting plasma are slightly larger 

than unity, suggesting that a beam is a poten¬ 
tially useful investment in these systems. 

The equilibration model employed here alters 
these conclusions somewhat. At low temper¬ 
atures, the pressure is reduced and the beam 
reactivity benefits from the reduced thermal 
broadening of the fusion resonance. R conse¬ 
quently improves, but the collisionality of the 
plasma increases simultaneously. Faster relax¬ 
ation rates require the diversion of additional 
fusion power to maintain the beam and reduce 
the possible temperature separation between 
species, increasing the effective bremsstrahlung 
emission. Both of these effects limit the eco¬ 
nomic viability of a fast beam. However, the 
model excludes beam particles from the calcu¬ 
lation of collision rates, reducing the collisional¬ 
ity of the plasma as the beam fraction is further 
increased. These competing collisional effects 
constitute the primary interactions of the beam 
with the thermal plasma. 

Apart from the beam fraction, the other 
salient parameter in scanning plasma conditions 
is the number ratio e, defined in all cases as the 
ratio of heavy ions to light ions, including beam 
ions. An interesting divergence arose between 
the two reactions: in DT, e > 1 plasmas had 
lower ignition thresholds in both MCF and IGF. 
In p-^^B , e < 1 plasmas were favored for MCF 
and e > 1 plasmas were favored for IGF. In DT, 
the number of electrons does not depend on e, so 
the bremsstrahlung emission is mostly indepen¬ 
dent of this parameter. The slight preference for 
tritium-rich plasmas is likely seen in extremal 
cases because the same center of mass energy is 
available at lower pressures in a deuterium-rich 
plasma. Inp-^^B , the preference for boron-poor 
plasmas in MCF is explained by the balancing 
of bremsstrahlung emission (increasing with e) 
and thermal reactivity (decreasing with e). 

The large boron density in the IGF optimum 
is harder to explain. The number ratio which 
maximizes the reactivity of a pressure-limited. 


constant-temperature system is e = = 

(1 + Zb)/(1 -I- Za) = 1/3, in the case of p- 
^^B. However, if the species temperatures are 
allowed to float, a high electron density (due to 
the boron excess) could suppress the electron 
temperature enough via bremsstrahlung emis¬ 
sion that the reduction in the thermal content 
of the plasma outpaces the loss of reactivity and 
improves the volume gain in the aggregate IGF 
accounting. 

In the presence of a beam, the e and Ti 
maximizing the total reactivity-pressure ratio 
should increase and decrease, respectively, as 
these changes create a rich environment of cold, 
resonant target ions for the lighter beam ions. 
The result is increased reactivity at fixed pres¬ 
sure. In p-^^B , the price of increasing e 
is rednced thermal reactivity and, moreover, 
greater bremsstrahlung emission. In a constant- 
temperature model, both fuels behave as pre¬ 
dicted; as the beam fraction rises from 0 to 1, 
there is a critical point where the beam reactiv¬ 
ity exceeds the total thermal reactivity and the 
optimum operating point is an excess of very 
cold target ions. 

Although the volume gain calculation is in¬ 
trinsically conservative (cf. ignition in a hot 
spot), the assembly energies Ea have been sys¬ 
tematically underestimated, thereby exaggerat¬ 
ing the gains reported here. In practice, it 
will be necessary to prepare the target at the 
prescribed temperature, which will likely in¬ 
volve heating from low temperatures (perhaps 
cryogenic, ambient at best). This heating will 
have to overcome bremsstrahlung and conduc¬ 
tion losses; a better estimate of the volume 
gain will require tracing the target’s evolution 
in pR — T phase space as Lindl has done for DT 
hot spots in NIF targets.!^ 

Investing substantial recirculated power in a 
fast beam has the potential to reduce ignition 
threshold conditions and ameliorate various en¬ 
gineering difficulties. The encouraging results 
of the stability analysis suggest that because 
the bulk ions are usually sufficiently cold, the 
necessary recirculated power can be regarded as 
small. In the case of DT plasmas, fast beams 
can increase the reactivity of a magnetically 
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confined plasma at constant pressure or lower 
the pR required for ignition by an inertially con¬ 
fined plasma. In plasmas, the extra elec¬ 

tron density increases the cost of maintaining a 
beam, limiting the gains which are possible in 
principle. 

Indeed, this work found superior gains in DT 
with the introduction of a beam. In both MCF 
and IGF DT plasmas, these gains were substan¬ 
tial. A fast deuterium beam improved the DT 
reactivity by an order of magnitude while reduc¬ 
ing the ion temperature to order-1 keV levels. 
However, the 10"^ reduction in the IGF igniting 
assembly energy is arguably even more impres¬ 
sive. In IGF plasmas, the correspond¬ 

ing reduction was only a factor of 30. How¬ 
ever, this reduction is more impressive than the 
gains realized in p-^^B MGF, where a beam im¬ 
proved the total reactivity by only 20%. How¬ 
ever, the alpha channeling which could be used 
to power the beam is likely critical to any steady 
state, magnetically confined scheme as a reliable 
means of abating the pressure poisoning effect. 

A key caveat to all the conclusions reached 
here, and the comparisons made between IGF 
and MFE as well as between DT and p-^^B fu¬ 
sion, is that these conclusions and comparisons 
are all based on an upper bound to a utility that 
in practice may be difficult to reach. However, 
the utility that can possibly be reached serves 
as impetus to try to find ways to reach it. The 
scenarios that we considered are in that respect 
at least not disallowed by the laws of physics, 
so they can serve at once not only as an impe¬ 
tus to achieve what is not disallowed, but also 
as a caution not to expect that more could be 
achieved. 
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Appendix A: Maintaining a monoenergetic beam 

There are several possible model choices for 
the recirculated power required to maintain the 
proton beam. RideJ^ calculated the amount 
of power necessary to prevent velocity space 
spreading of the beam past a prescribed thermal 
width vth ^ Vb- This spread results primarily 
from beam-beam collisions; the Rider estimate 
of the recirculated power is proportional to n^. 
This calculation is applicable when the beam 
lies far outside the thermal proton population 
and the velocity width of the beam is impor¬ 
tant. 

In general, Precirc must also counteract beam 
collisions on the thermal population, involving 
terms proportional to ribUp, ngriB, and nyrie. 
A more precise calculation would include both 
beam-beam and beam-thermal collisions and 
evaluate the reactivity integrals with a realistic 
slowing down distribution for the beam protons. 
In pursuing an upper bound of the gains realiz¬ 
able from alpha channeling, it has been assumed 
that the monoenergetic beam collides only with 
thermal particles. 

1. Collisional evolution of beam 

The Fokker-Planck equation governs the col¬ 
lisional evolution of an initially monoenergetic 
beam. Drag and velocity space diffusion in¬ 
duced by each of the thermal species contribute 
to the beam’s eventual thermalization. We in¬ 
voke a test particle analysis to argue for the 
validity of the model employed here. Consider, 
for example, the case of a 600 keV proton slow¬ 
ing in a plasma with Tp = Tb = 120 keV and 
Te = 80 keV. The number ratio is ub/pp = 0.20 
and rie = 3.5 x I0^^cm“^. Although the beam 


slows primarily on thermal protons until an en¬ 
ergy scattering time has passed, the fusion rate 
falls significantly from its original value on a 
collisionless time scale. 

A similar DT case would have a 125 keV 
deuteron slowing in a plasma with To = Tt = 
llkeV and = 24keV. The number ratio is 
pt/pd = 1.25 and Pe = 7.5 x 10^^cm“^. The 
picture is altered here because the ideal beam 
energy (maximizing the beam-thermal reactiv- 
is somewhat larger than the fusion reso¬ 
nance, such that most of the beam fusion events 
come after about one collision time. However, 
the fusion rate remains larger than the thermal 
ion collision times throughout the slowing down 
process. 

In all cases, the beam remains well outside 
the thermal ion distributions, and most fusion 
events occur on a collisionless timescale. We do 
not anticipate that collisional spreading of the 
beam will modify the power balance laid out 
here. 


2. Collisionless instability 

Let us expand on comments made in the 
Introduction concerning the stability of the 
monoenergetic beams. In particular, we will 
present suggestive criteria for ion beams propa¬ 
gating in either cold or warm plasmas. A more 
detailed stability analysis for a collisional three- 
component plasma is beyond the scope of this 
work. 

In a collisionless two-component plasma, the 
stability threshold for Langmuir waves in the 
presence of an ion beam igiTI 

Vlu) = Z'(u)po^Z'[a{u — V)] (Al) 
-2{p + l)T = 0, 

where Z' is the first derivative of the plasma dis¬ 
persion functiorP^, u = ui/kvth,i is a dimension¬ 
less phase velocity, a = Vth,i/vth,b describes the 
‘coldness’ of the beam, V = Vblvth,i measures 
the beam’s separation in velocity from the bulk 
ions, and p = Pbjni is the dimensionless beam 
density, and T = Ti/Tf.. Generally speaking. 
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V is Oil) in all configurations and T is larger 
forp-^^B than DT plasmas (greater electron-ion 
temperature separation). 

Taking asymptotes of Eq. (Al), we can get 
some general notions of the stability of the 
beam-plasma conhgurations to be presented.^ 
We will present the stability threshold as a 
curve in (V,a) space parametric in the phase 
velocity u. For V ^ 1, corresponding to a 
beam well outside of the thermal ion bulk (cold 
plasma), the stability threshold is 

and likewise for E <C 1, a beam embedded in 
the thermal bulk (warm plasma), 

V = + (A3) 


The minimum thermal spread is then 1 keV, 
compared to a beam energy of about 627 keV. 

In these representative examples, quasi- 
monoenergetic beams are well within the sta¬ 
bility threshold. However, as the beam fraction 
is increased, the necessary stabilizing thermal 
spread may approach the beam energy in some 
scenarios. 


Appendix B: Optimized reactivity subject to 
maximum pressure 

In pressure-limited systems, the greatest 
achievable reactivity is an important experi¬ 
mental parameter that establishes confinement 
criteria. In the absence of a fast beam, the max¬ 
imum of the function 


with 


In practice, these criteria may be regarded as 
a minimum amount of thermal spread intrinsic 
to the beam. Cold beams with little thermal 
spread (large a) lying above these curves are 
unstable. 

For example, we consider a case typical of 
the configurations to follow, a DT plasma with 
n = 0.1 and T = 0.5. In both DT and 
plasmas, F > 1, so we will proceed cautiously 
using the F S> 1 asymptote. Then a stable 
beam has 

a<VnV-l. (A5) 

For a representative F = 4, a < 4-\/TT— 1 « 12, 
corresponding to a minimum beam temperature 
for stability of 100 eV in DT. (Typical beam 
energies are in the neighborhood of 127 keV, 
the center of mass energy at the fusion cross 
section peak.) Compare a typical case in p- 
, with n = 0.1 and T = 2. We assume 
the beam-plasma interaction largely ignores the 
boron ions. The corresponding stability crite¬ 
rion is a < VllOF — 1 « 20 for a typical F = 2. 


UA-riB ((Tv){T)\nAT+ nBT+ neT <po (Bl) 


is of interest, ua and ub are the number 
densities of the reacting ion species, such that 
Ue = Zatla + ZbUb- Po denotes the maxi¬ 
mum allowed pressure and UAfnB = e. Equiva¬ 
lently, one can maximize the Lagrange function 
L = UAnB (crv) — X{nAT + ubT — pq). The 
components of the gradient are 


- (1 + Zb + e(l + Za))TX + 2^^^^ 

(B2) 

— jr = iZB- Za)TX+ -——^rie {av} 

(B3) 

77-e ul H” ^ 


dL 1 -|- Zb + e(l + Za) „ 

AT = -T—“ Po 

UA Zb + Zj\€. 


(B4) 

(B5) 


Considering the complementarity condition, 
A = 0 implies net {crv) = 0, which cannot hold 
in a finite-temperature plasma. Proceed, fixing 
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p = Po- Then 


He 


^ 1 + Zb 
1 + Za 

(Za + Zb + 2ZaZb) Po 


2{1 + Za){1 + Zb) T 

i.vym = 2MT> 


where n = n(t) is the cumulative number of 
binary fusion reactions. Defining (p = 2n/no, 


(B6) 

dn 

(B7) 

dt 

no 

(B8) 

4(1+ e): 


= y-^ =(crv) (C5) 
((e + 1 ) 4 ) — 2)(e((/> — 2) + ^) 


where the final equation can be solved implicitly 
for the T = Tq which maximizes the reactivity. 

Appendix C: Evaluation of bum fraction 4> 


which is a Riccati equation for 4>- Assuming 
constant Tg and Ti during burn, integrate over 
the confinement time, R/iCg, to obtain the im¬ 
plicit equation 


Without loss of generality, suppose initial 
number densities uao and ubq such that n^o + 
i^Bo = no and nAo/nBo = £■ Then 

nAo = (Cl) 

1 -b e 

ubo = ^ ^ no (C2) 


1 + e 


1 - e 


log 


e(^(e -I- 1) — 2) 

e{4> — 2) -b (/) 


■:^no {av) 


(C6) 


This can be inverted to obtain 


and 


HAit) = p^no - n 
1 -b e 

nB{t) = - n 


(C3) 

(C4) 



1 - 


e- 1 


(C7) 


where p = niAnAo + wsnBo, Hb = 6csm/{av), 
and m = (eruA + nriBpli^ + 1). 



